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In this paper, we study Gibbs point processes involving a hardcore interaction which is not 
necessarily hereditary. We first extend the famous Campbell equilibrium equation, initially pro- 
posed by Nguyen and Zessin [Math. Nachr. 88 (1979) 105-115], to the non-hereditary setting 
and consequently introduce the new concept of removable points. A modified version of the 
pseudo-likelihood estimator is then proposed, which involves these removable points. We con- 
sider the following two-step estimation procedure: first estimate the hardcore parameter, then 
estimate the smooth interaction parameter by pseudo-likelihood, where the hardcore parameter 
estimator is plugged in. We prove the consistency of this procedure in both the hereditary and 
non-hereditary settings. 
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1. Introduction 

Gibbs point processes first appeared in statistical physics for the description of large 
interacting particle systems. They are now used in many other fields such as biology, 
medicine, agronomy and astronomy. Gibbs models rely on interaction potential functions 
that must be properly estimated in applications. We are particularly interested in the case 
where the potential function involves a hardcore interaction. This means that some point 
configurations are forbidden with respect to this potential function or, equivalently, that 
some point configurations have infinite energy. This constraint appears in many classical 
models of Gibbs measures (for example, the hard ball model). 

To the best of our knowledge, in all previous studies involving a hardcore interaction, 
the interaction is assumed to be hereditary (see [15]). As recalled in Definition 1, an 
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interaction is hereditary if, for every forbidden point configuration 7 in R and every x 
in R'^, the configuration j + 6x (that is, the configuration composed of the union of the 
point X with the points of 7) remains forbidden. Note that any interaction involving no 
hardcore part is necessarily hereditary. 

In this paper, we aim to study Gibbs measures in the presence of a hardcore in- 
teraction that is not necessarily hereditary. Indeed, it seems natural to encounter some 
non-hereditary interactions, as Examples 5.2 and 5.3 in Section 5 show. The first example 
is concerned with random fields of geometric objects subjected to a hardcore interaction. 
Such a model is developed in [5] with the study of random tessellations with geometric 
hardcore interaction. In this model, each cell of the tessellation is forced to have a radius 
less than a fixed constant a. It is clear that this kind of interaction is not hereditary. In 
fact, when one adds a new point inside a forbidden large cell, the new tessellation may 
become authorized. The second example is a forced clustering process. In this model, 
each point is forced to have at least k neighbors at a distance less than a. Here, again, 
this interaction is clearly not hereditary. 

In the first part of the paper, we introduce some notation, recall the definition of 
Gibbs measures on M.'^ and prove some preliminary results about non-hereditary Gibbs 
processes. The main problem, in the non-hereditary case, is that the energy of a point x 
in a configuration 7 is not always defined. Indeed, the energy of "f — dx may be locally 
infinite, even if the energy of 7 is locally finite. In this case, the energy of a; in 7 would 
be negative infinity, which makes no sense. We must then introduce the set TZ{'y) of 
removable points in 7 (see Definition 2): x in 7 is said to be removable if the energy of 
"f ^ Sx is locally finite. 

The second part of the paper is dedicated to the elaboration of a Campbell equilib- 
rium equation for non-hereditary Gibbs measures. The classical Campbell equilibrium 
equation, proposed by Nguyen and Zessin [17], is only valid in the hereditary case. This 
equation is fundamental for the understanding of Gibbs measures. It also provides major 
statistical applications: the Takacs-Fiksel and pseudo-likelihood estimation procedures 
rely on the Nguyen-Zessin equation (see, for instance, [15] for an exposition). Moreover, 
the recent definition of residuals for spatial point processes arises from this equation (see 
[3]). The Introduction of the concept of removable points allows us to extend the Camp- 
bell equilibrium equation when the heredity of the interaction is not assumed. Let fi be 
a Gibbs measure and C'^ its reduced Campbell measure. We prove in Proposition 2 that 

(d^- d7) = e-''^-''^) A ® fi{dx, dj), (1) 

where h denotes the local energy of 2; in 7 and A is the Lebesgue measure on W^. 

Next, we consider the estimation problem. Several statistical procedures have been 
proposed to estimate the potential function of a hereditary Gibbs measure. Among them, 
the common parametric methods are the maximum likelihood estimation and the pseudo- 
likelihood procedure. In addition, the Takacs-Fiksel method, some Baycsian procedures, 
and semi- and non-parametric methods have been developed. We refer to [15] for a recent 
review of these estimation procedures. 

In Section 4, assuming the interaction is parametric, we focus on the pseudo-likelihood 
method. The direct parametric method consists of the maximum likelihood procedure. 
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but it requires the estimation, by simulations, of an unknown normalizing constant. 
Recent advances, including Monte Carlo and perfect simulation, deal with these simu- 
lations efficiently (see [4, 7, 13]). Yet, from a practical point of view, maximum likeli- 
hood remains computationally much more intensive than pseudo-likelihood. Therefore, 
pseudo-likelihood estimation is often used as a first step before maximum likelihood es- 
timation. For the same reason, it may be preferred over maximum likelihood estimation 
when fitting several alternative models to a data set. Moreover, in the stationary case, 
some asymptotic results exist for the maximum likelihood procedure (see the review in 
[15]), but they rely on rather restrictive assumptions that seem difficult to avoid. The 
pseudo-likelihood method has the advantage of satisfying the usual asymptotic theory 
for a large class of interactions; see [2, 8-10, 14]. In these articles, under very general 
assumptions, the consistency of the pseudo-likelihood method is proved, as well as the 
asymptotic normality of the induced estimators. In the same spirit, in Section 4, we prove 
the consistency of the pseudo-likelihood procedure when the interaction is not necessarily 
hereditary. 

More precisely, we suppose that the interaction depends on two parameters, a and 
0. The parameter a is devoted to the hardcore interaction, whereas 6 parameterizes the 
classical smooth interaction. In our non-hereditary setting, the pseudo-likelihood contrast 
function is defined as 



|A„| 



(2) 



where A„ denotes the observation domain of the sample. This pseudo-likelihood contrast 
function is the classical one found in previous studies, but in order to remain well defined 
in the non-hereditary setting, the sum is restricted to the removable points 7^" (7) (this 
set depends only on the hardcore parameter a, as implied by assumption SI in Section 4). 
We consider a classical two-step estimation procedure. We first estimate a in a consistent 
way (see Proposition 3). We then estimate 9 by maximizing PLL\^, where a is replaced 
by its estimator. Theorem 2 establishes the consistency of both estimators of a and 9 
resulting from this procedure. Let us note that even if the interaction is hereditary, to 
the best of our knowledge, no consistency proof has previously been supplied for this 
classical two-step procedure, although it has been widely used in practice. This is, for 
instance, the case for the classical hardball model that we present in Section 5. 

Section 4 is devoted to the consistency of the estimators, but their asymptotic normal- 
ity is not addressed. However, there are no major difficulties in obtaining the asymptotic 
normality of the estimator of 9 when a is known, provided the interactions are finite 
range. The tools for proving this are essentially the same as in [10] and, more recently, 
in [2]. In the situation where the hardcore parameter a is estimated, the asymptotic 
normality of the estimator of 9 seems more difficult to obtain. 

Section 5 contains some examples of parametric interactions involving a hardcore part. 
The first is the hardball model with pairwise step interaction. It is a hereditary model. 
The second is the hardcore Delaunay tessellations model, which is not hereditary. The 
third example is a forced clustering process with interaction, which is also not hereditary. 
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For these examples, we prove that the assumptions of Theorem 2 are satisfied, leading 
to the consistency of the joint estimation of their parameters. 

2. Definitions and notation 

2.1. State spaces and reference measures 

Let d be a fixed integer greater than 1. Mf^ denotes the d-dimcnsional Euclidean real 
space endowed with the Borel cr-algebra cr(]R'^). B(M'^) is the set of bounded Borel sets 
on R''. y\4(R'^) is the set of simple integer-valued measures 7 on M.'^, that is, for every 
A e B(R''), 7(A) G N and for every x G R'^, 7({a;}) < 1. X(R'*) is endowed with the 
cr-algcbra ct(X(R'')) generated by the sets {7 G 7W(R''),7(A) = n}, n G N, A G S(R''). 
So, V(M{M.''^)) denotes the set of probability measures on AdiM!^) for the cr-algebra 
a{A4.{W^)). Any measure 7 G Ai{Mf'-) has the following representation: 

where I is a subset of N, {xi)i^x are elements of R'' and Sx is the Dirac measure at x. 
We write a; G 7 if 7({a;}) > 0. 

Letting 7 be in A^(R'^) and A be a Borel set in R'^, we denote by 7a the projection of 
7 on A, which is just the measure J^x^-ynA^^- 

Denoting by A the Lebesgue measure on R'', tt stands for the Poisson process on R** 
with intensity A. It is a probability measure on A^(R''). For every A G i3(R'')), tta denotes 
the Poisson process on A with intensity 1a A. It is a probability measure on M{A). 

2.2. Interaction 

We define the interaction energy in a general setting, as in [12]. We consider the most 
general specifications of Gibbs kernels with local densities having an exponential form. 
We do not assume that the local densities come from a multibody interaction potential. 
This general point of view allows us to deal with the non-hereditary case (sec Definition 
1 below). 

A set of measurable functions (^^A)AeB(R'') from A^(R'') to KU {-|-cxd} defines a family 
of energies if, for every A C A' in i3(R'^), there exists a measurable function (^a,A' from 
M{M.'^) to RU {+00} such that 

V7G>I(R'^) i/A'(7) = ffA(7) + '^A,A'(7A0- (3) 

(3) is equivalent to (6.11) and (6.12) from page 92 of [12]. In particular, we have the 
following property: 

VAcA' in6(R'') andV7GX(M'*) Ha{"/) ^ +00 => i7A'(7) = +oo. (4) 
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Physically, H\{'-^) is the energy of 7a inside A given the configuration outside A. Let 
us now discuss the problem concerning heredity. 

Definition 1 . A family of energies is hereditary if 

VAeS(M''),V7eX(R'') andVa-e A Fa(7) = +oo ^ Ha{i + 5x) = +oo. (5) 

This assumption (5) is often needed in many papers (see, for example, [16, 17]). Let us 
point out that the non-heredity of a family of energies necessarily comes from a hardcore 
part (that is, configurations with infinite energy). From now on, we will never assume 
that the energy is hereditary. When we invoke the "non-hereditary case" , it means the 
general setting of hereditary or non-hereditary energies. 

Let us introduce the new concept of removable points, which is fundamental in the 
non-hereditary case. 

Definition 2. Let 7 be in M{W^) and X he a point of"/, x is removable from 7 if 

BAeBiR'') such that xe A and Ha (7 - 4) < +00. (6) 
We denote by K(7) the set of removable points in 7. 

When the configuration 7 has a locally finite energy, this definition can be simplified, 
as proved in Proposition 1. A configuration 7 in M-(W^) has a locally finite energy if 
for every A e B(M.'^), the energy H\{j) is finite. We denote by Mooi^'^) the space of 
configurations which have a locally finite energy. 

Proposition 1. If -f is in A^oo(II^^) o,nd x is a point of j, then x is removable from 7 
if and only if J — 6x is in A4oo(^^)- 

Proof. Let 7 be in A^oo(I^'*) and a; be a removable point of 7. There exists A e S(M'*) 
such that a; is in A and H^^j — Sx) is finite. Let us show that -ffA'(7 — S^) < +00 for 
every bounded set A' in R'^. Thanks to (3), we have 

Ha'ua{1 - Sx) = HAi"/ - Sx) + </'A.AUA'((7 " Sx)a'') 
= Ha{-1 - Sx) + 'PA,AUA' ilA- ) ■ 

Since 7 is in Aiaoi^'^), 'Pa,auA' {jA") is finite. Thus, -ffA'uA(7 — Sx) is finite and, thanks 
to (4), Ha' (7 - Sx) is also finite. Therefore, j — Sx is in A^oo(R'^)- 

The converse is obvious. □ 

We can now define the local energy of a removable point x in a configuration 7. 

Definition 3. Let x be a removable point in a configuration 7 in y\4(R'^). The local 
energy of x in "j — Sx is defined as 

h{x,"i-Sx)^HA{i)-HA{i-Sx), (7) 
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where A is a bounded set containing x such that H\{'-^ — Sx) is finite. 

Notice that this definition is vahd and does not depend on the choice of A. Indeed, 
according to Definition 2, at least one such A exists. Besides, let us suppose that there 
is another A' containing x and such that H\r{'y — 5x) is finite. Writing A" = A n A', we 
have, from (3), 

i?A(7) -Ha{j-S.j;) = i/A"(7) + <PA",a(7A"c) - HA>'ij - Sx) - VA",a((7 - '5x)a"0 
= i7A"(7) + 'i5A",a(7A"0 - HA"ij - 5x) - <y5A",A(7A"0 

= Fa"(7)-^A"(7-^x), 

which is equal, thanks to the same calculations, to i?A'(7) ~ ^^A'(7 ~ ^x)- Therefore, 
h{x, 7 — Sx) is well defined and belongs to R U {oo}. 

Finally, for every 7 in A4oo(R'') and every x in R'', the local energy of a; in 7 is always 
defined since h{x^ ^) — h{x^ ^ '\' ~ 

) and Ha (7) is finite for all bounded Borel sets A 

in W^. 



2.3. Gibbs states 

In this subsection, we are in a position to define the Gibbs states. We need to introduce 
the notion of specifications, as in [12]. Let us make an integrability assumption about the 
family of energies, which is equivalent to (6.8) in [12]. We say that the family of energies 
{Ha) is integrahle if, for every A in S(R'^) and every 7 in 7Woo(R''), we have 

0< / e-^^('^A+'^^'^)7rA(d7;v) <+oo- (8) 

Jm{K) 

The second inequality in (8) is, in general, guaranteed by the stability of the potential. 
We will assume this stability in H3 below. The first inequality is obvious in the classi- 
cal hereditary setting. In the non-hereditary setting, it remains true under reasonable 
assumptions (see, for instance, [5]). 

Under this integrability assumption, we are able to define the kernels for the Gibbs 
structure (see (6.6) in [12]). For every A in B(R''), the kernel Sa on 'P{Moo{^'^)) x 
Moo{^'^) is defined by 

Sa(/,7) = / /(7')SA(d7',7) 

(9) 

:= / /(7; + 7Ac)— ^e-^'^(^A+^-)7rA(d7;), 

JM{A) -^A(7A-j 

where / is just a bounded measurable test function and Za{'^A':) is the normalization 
constant defined by 

Za{im)= [ e-^^(''A+''A=)7rA(d7;). 
Jm{a) 
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Note that from (8), < Z\{'y\c) < +oo and therefore the kernels are well defined. More- 
over, due to (3), they are compatible in the sense of (2.14) in [12], which means that for 
every A C A', every 7 in Moo(M'') and every bounded measurable function /, 

/ /(7')SA'(d7',7)= / / /(7")5A(d7",7')5A'(d7',7). (10) 
Now let us give the definition of Gibbs measures (see (2.15) in [12]). 

Definition 4- A probability measure ^ on A4{R'^) is a Gibbs measure for the family of 
integrable energies (Ha) if, for every A in B{M.'^) and every bounded measurable function 
f from MiW^) to R, we have 

I fil)Kdl)= f f /(7')5A(d7',7)M(d7)- (H) 

We denote by Q the set of stationary Gibbs measures. Equations of the form (11) are 
called DLR equations, where DLR stands for Dobrushin, Landford and Ruelle. They may 
be rewritten in the following way: for /x-a.e. 7 and for every bounded set A in B{W^), 

/i(-|7Ac) = Sa(-,7). 

Remark 1. From (11), we deduce that the support of /i is included in A^oo(IR'^)- 

3. Campbell equilibrium equation for non-hereditary 
Gibbs point processes 

In this section, we develop a Campbell equilibrium equation for non-hereditary Gibbs 
point processes. In [17], the authors give an equation using the reduced Campbell measure 
to characterize hereditary Gibbs point processes. This Nguyen-Zessin equation (12) is 
very well known and is used in many works concerning Gibbs processes: 

C|,=c-''A®M- (12) 

Unfortunately, this formula is not valid in the non-hereditary case, as explained in 
Remark 2 below. Therefore, we propose to generalize it, in Proposition 2, to the non- 
hereditary case. In the hereditary setting, it simply becomes the classical Nguyen-Zessin 
equation (12). The concept of removable points introduced earlier is the key for this 
generalization. 

First, let us introduce the definition of the reduced Campbell measure (see, for example, 
[11], page 225). 
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Definition 5. Let n be a probability measure on MiR''). The reduced Campbell measure 
Cl on M'^ X A^(K'^) is defined by 

C-(/)= / / /(x,7-4)7(dx)M(cl7), 

where f is a bounded non-negative measurable function from x A^(]R'*) to R. 

Remark 2. We shall explain why the Nguyen-Zessin equation is not valid in general for 
a non-hereditary Gibbs point process. Let us show that is not absolutely continuous 
with respect to A(8)/i, which implies that (12) is not satisfied. For this, since the support 
of ji is A^oo(IR'^), it is sufficient to prove that the support of is not included in 
M'' X Mao{W^). Indeed, from Definition 5, 

CliiW'xMU^'W)^ I I V*-.^A^oo(K^)7(dx)M(d7), 

JM(R'^) JRi^ 

where \A\'^ denotes the complementary set of A in R'' x AA{W^). But, for a non-removable 
point a; in 7, 7 — is not in Moo{^^^) (see Proposition 1). Since non-hereditary Gibbs 
point processes contain some non-removable points, C^([R'^ x A^oo(R'')]'^) does not vanish 
in general and the support of is not included in R'' x M.ao{W^)- This is exactly the 
situation in Examples 2 and 3 in the last section. 

Let us now present our generalization for the non- hereditary case. 

Proposition 2. Let ^ be a Gibbs measure in Q . For every bounded non-negative mea- 
surable function f from x A^(R'') to R, we have 

I lAioo(R'')(7)/(a;,7)C),(dx,d7) 

JR<^y.M(R<^) 

(13) 

= / / /(x,7)e-''(-'^)A(dx)M(d7). 

JR-^ JMiR'^) 

We have seen in Remark 2 above that, in general, the support of C'^ is not included 
in R'^ X 

A4oc(^'')- The Campbell equilibrium equation (13) shows that restricted to 
M!^ X A^oo(K'') is absolutely continuous with respect to A®^ with density e^''. 

Let us point out that this proposition is also valid in the non-stationary case. The 
intensity A is then replaced by any locally finite measure. 

Moreover, it is important to note that the converse of Proposition 2 is not true, which 
means that (13) docs not characterize the measure /i. Consider, for example, a measure 
fi such that, almost surely, 7 docs not contain any removable points. (13) then becomes 
the trivial equation = 0. In fact, the equilibrium equation (13) is interesting only if, ^ 
almost surely, 7 contains some removable points. 
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Proof of Proposition 2. Let /i be a Gibbs measure and / a bounded non-negative 
measurable function from R'^ x A^(M'') to M. Let A be a bounded set in M'*. By tlie 
definition (5) of tlie reduced Campbell measure, we have 

/ lA^<^(K'')(7)lA(a;)/(a;,7)C|,(da;,d7) 

Jm{R''} Jr'' 

By the definition of the local energy, Proposition 1, equation (9) and the DLR equations 
(11). we have 

/ ^M^(R-^)il)'^Aix)fix,j)Cl{dx,d-f) 

jR'^xMiR'') 

= V lA(a;)/(a;,7-4)^(d7) 

•^^(«'').eK(7)nA 

= [ [ E /(^,7'-'5.)SA(dy,7)M(d7) 

= / / E /(^.7;+7A--'^.) ^( . 7rA(d7;Md7) 

•^^(R'')-^^(A),eK(7;+7Ao)nA ^^^^^^^ 

= 11 V /(x,7;+7Ac-4)e-''(^^^^+''^-'-) 

JA1(R^)JA1(A)^g^(^,^^^^)^^ 

g-ifA(7A+7A<:-l5x) 

^ Zd^) -A(d7;)Md7) 

= / / El^(7U7Ac+5.-^.)(^)/(2;.7U7Ae-<5,)e-''(^'''A+^--*^) 

g--tfA(7A+7Ac-5x) 

— ^ira — "^^^^^^'^^'^^^ 

= / / l7^(^;+^A.+5.)(x)/(x,7;+7A0e-''^^^^^+^-) 

g--HA(7A+7A=) I 

^a(7aO ^^A(d^'d^A)M(d7). 
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The well-known Campbell equilibrium equation for the unit rate Poisson process (that 
is, Cl^= X^itta) gives 

[ lniY,+,,.+s.){x)f{x,j'A+iM)e-^^^'^'^+''^^^ 

A1(K<*) JAxM(A) 

g--H"A(7A+7Ac) 

X — 7n ^ — A(g)7rA(dx,d7A)//(d7). 

Again using (9), we have 

/ lAi<»(R<i)(7)lA(a;)/(x,7)C|,(dx,d7) 

= / / / lK(y+,^)nA(x)/(x,7')e-''(--^')A(dx)SA(d7',7)A^(d7)- 

JR'^ JMiR'') JMiR'') 

Note that if 7' is in A^oo(K'') and x is in M.'^, then x is in TZ{j' + Sx). Therefore, we have 
/ ^M^(R-^)h)'^A{x)fix,-f)Cl{dx,d'y) 

Jr'^xM(R'') 

= 11 I lA(x)/(x,7')e"''^"''''^A(dx)SA(d7',7)M(d7) 

jR'i JM(R'i) JMiR-^) 

= 11 lA(x)/(a;,7)e-^(^''^U(dx)Ai(d7). 

jRrf JMiR-^) 

Relation (13) is proved for every A in S(M^), so it is proved for A = M'' as well. □ 

4. Consistency of the pseudo-likelihood estimator 

In this section, we deal with the parametric estimation of (non-hereditary) stationary 
Gibbs measures. Our aim is to prove the asymptotic consistency of the estimation pro- 
cedure when the observation window of the Gibbs point process increases to R'^. We 
focus on the pseudo-likelihood procedure, which appears to be asymptotically validated 
for a large class of interactions. This procedure, quicker in practice than the maximum 
likelihood approach, may constitute a first consistent estimation in applications and may 
therefore be used to quickly fit several alternative models. 

We suppose that the family of energies {H\) depends on a positive parameter 0* and on 
a multiple parameter 6* = {6\, . . . ,9*^). The first parameterizes the support of the energy 
(that is, when the energy is equal to positive infinity), while the second parameterizes 
the energy when it is finite. These two parameters play very different roles, as the set of 
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assumptions made in this section will make clear. Note that it would be easy to consider 
a vectorial hardcore parameter, but, for reasons of clarity, we chose to focus solely on a 
real hardcore parameter. 

Let {a, 9) be two parameters in K+ x 9, where 6 is a bounded open set in W. Let us 
denote by (J^a'^)agB(r<^) the parametric family of energies and Q"'^ the set of stationary 
Gibbs measures for this family of energies. The following assumption SI implies that the 
support of the energy is parameterized by a only, not by 9. 

SI. For all 7 £ 7W(M'^), a e R+ and 9, 9' in 9, 

^keB{W) Hl'\-f) <oo ^ Hi'''' {-f)< 00. 

This assumption claims that the set of configurations 7 in 7W(M'') which have a locally 
finite energy for the family (i?^' ) depends only on a. We therefore denote this set by 
A^^(IR''). The same remark is true for the set of removable points in 7, and we similarly 
denote this set by R"(7). Finally, for every x in R"(7), we define /i"'^(x,7 — S^)-, the 
energy oi x hi ^ — 5x, as in (7). 

For all a G and 6' G 9, wc define the pseudo-likelihood function at 7 G M.'^iW^) as 



PLLA„(7,a,^) = -i 



cxp(-/i"'''(x,7))da;+ ^ h"''ix,j-Sx) 



(14) 



where A„ denotes the domain of observation of the sample and |A,i| is its Lebesgue 
measure. 

The purpose of this section is twofold. First, when the hardcore parameter is supposed 
to be known, we aim to prove the consistency of the pseudo- likelihood estimation of 9* , 
without assuming the heredity of the interaction. This problem is addressed in Section 
4.1. Since we are in a non-hereditary setting, the definition (14) of the pseudo-likelihood 
contrast function differs from the common one: it involves the set of removable points. In 
the hereditary case, the proof of the consistency relies mainly on the ergodic theorem and 
the Nguyen-Zessin equation. In the non-hereditary case, since the Campbell equilibrium 
equation is modified into (13), we must check carefully the consistency of the pseudo- 
likelihood procedure. We prove in Theorem 1 that the non-hereditary setting does not 
actually modify the range of validity of the consistency since the set of hypotheses on 
the smooth interaction parameter 9 is the same as in the hereditary case (for example, 
as in [2]). 

The second purpose, addressed in Section 4.2, is to prove the consistency of the two- 
step estimation procedure, without assuming the heredity of the interaction. The support 
parameter a* is first estimated and plugged into the pseudo-likelihood contrast function. 
9* is then estimated by the common pseudo-likelihood procedure. To the best of our 
knowledge, this two-step procedure, although widely used in practice, has not yet been 
proven to be consistent, even in the hereditary setting. Apart from the usual conditions 
on the smooth interaction parameter, some additional assumptions on the support pa- 
rameter a are needed (see S2-S4): roughly speaking, they suppose that the sets A1J^(R'') 
are embedded and continuous with respect to a. 



Non-hereditary point processes 



1379 



In Section 5, some examples of models involving both a hardcore and a smooth inter- 
action are presented. We prove that they satisfy all the assumptions of this section. 

4.1. Consistency of 6^ when the support parameter a* is known 

Let us suppose that the first parameter a* is known. Wc estimate 6* by 

§n = argminPLLA„ (7, a*, 6). (15) 

eee 

To prove the consistency of this estimator, wc need the following set of hypotheses: 

HI. (A„)„>i is an increasing sequence of convex and compact sets such that | A„| — R''; 
H2. the energy function h is invariant by translation, that is, for all x and y in M'' and 
all 7 G A4(M''), h{x + y,^y) = h{x,^), where 7j, is the configuration 7 translated by y; 
H3. 3K > such that for aU {a, 9), 7 e M'^^iR'^) and x E W^, 

H4. for all 9 in 6, 

h'^'''io,■)cM-h"''''{o,■))eL\^l'^''''), 

with the convention ooe~°° = 0; 
H5. for aU 9 me\{e*}, 

^"'■'''(/i"*'^'(0,-)?^/i"''*(0,-))>0; 

H6. one can find a real function 6 with S{x) when a; — and g S ), such 

that V(6i,6i') e 92, V7 e M^iR"^), if /i"'^''(0,7) < +00, then 

|/i"''«(0,7) - (0,7)1 < .9(7)^(1^ - ^'D, /i"*''*-a.e. 

HI is a natural assumption concerning the domain of observation. H2 and H3 state 
that h is invariant under translation and is locally stable. Local stability is a stronger 
assumption than the classical stability hypothesis encountered in statistical mechanics, 
yet it is fulfilled in many Gibbs models (see, for instance, [1, 5]). H4 is a technical 
integrability assumption which holds in most models (see [1, 5, 14]). H5 guarantees 
that is a proper parameter for the energy. From H5, we deduce that, ^" -a.s., there 
exists some configuration 7 such that h" (0, 7) < +00. One could then prove, using the 
ergodic theorem, that 7 contains almost surely some removable points. As a consequence, 
the sum involved in (14) is non-empty. Finally, H6 specifies the sense in which 9 1— >■ 
/i" '^(0,7) is continuous at 9*. 

Most of these hypotheses are similar to the assumptions found in [2] . Indeed, the latter 
deal with the pseudo-likelihood estimation in the general hereditary case. We adapt their 
scheme to the non-hereditary setting. The assumptions and the proofs are, therefore, in 
the same spirit. 
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Theorem 1. Let yu" € G" ■ Under SI and H1-H6, the estimator On defined by 

(15) is strongly consistent, that is, fi" -a.e., 

lim dn==0*. (16) 

n— f oo 

The pseudo-likelihood procedure is a minimum contrast estimation. This point of view 
has been used in Jensen and Kiinsch [9] and in Billiot, CoeurjoUy and Drouilhet [2] to 
prove the consistency and asymptotic normality of their estimator. It relies on Theorem 
3.4.3 of Guyon, established in [8]. 

First, note that we only have to prove Theorem 1 for ergodic measures /i" . If 
is not ergodic, it can be represented as a mixture of ergodic stationary Gibbs measures 
(see [12]). Therefore, from now on, is assumed to be ergodic. 

The following Lemmas 1 and 2 allow us to apply Theorem 3.4.3 in [8], which yields 

(16) . 
Let 

Kn{e, e* ) = PLLa,, (7, a*, 9)- PLLa„ (7, a*, 9*). (17) 
We prove below that Kn is a proper contrast function. 

Lemma 1. Under SI and H1-H5, for all 9 and /i" -a.e., 

lim K„{9,9*)=K{9,9*), 

n— f oo 

where K{-,9*) is a deterministic positive function which has a unique minimum at 9* . 

Proof. The main point of the proof consists of proving that for all 9 £ Q and /.t" - 
almost every 7, 



lim PLLAj-f,a*,9) 

= E^^^e- [exp(-/i"*^^(0, 7)) + h,,, ,e(0, 7) exp(-/i"* (0, 7))] 

where Ea*fi-' denotes the expectation under /.t" . 

Thanks to II1-H3, we can apply the ergodic theorem (see [17]): 



(18) 



exp(-/i" •^(x,7))d2; 

[0,1]'' 



lim / cyi]i{~h°-'^%x,-i))dx^Ea',e' [ 

n-yoo |A„| [Jfj 

The stationarity of fi" yields 

lim f exp(-/i"^^(x,7))dx = i;„..e4exp(-/i"*''(0,7))] 

n-i-oo |A„| Ja.„ 

which proves the first part of (18). 
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For the second part, note first that from Proposition 2 and the stationarity of /.i" , 

From Definition 5 of the reduced Campbell measure, we deduce that 



a;G7?,°*.<'(7)n[0,l]'' 



which is finite, thanks to H4. Therefore, we can again use the ergodic theorem and 
Ai"*'^'-a.e. 



iim -— - 

a:G7?.°*(7)nA„ 



■xeK°*(7)n[04]'' 



= C|^o-,e* (l7Hc^*(Md)(7)^a*,e(2^,7)lxe[0,l]'') 

The convergence (18) is proved. 

From the definition (17) of Kn, we deduce from (18) that fi" -a.s., 

Iim K„ie,e*) 

n— j-oo 

= i?„.,,.[e-''°*.''*(0'^)(e'^-.-(0'^)-''°*.«(°''') - [/ia.,e.(0,7) - h^.^O,!)] " 1)] 

According to H5 and the behavior of 1 1-)- e* — i — 1 , wc conclude that K is positive and 
vanishes if and only if 61 = 6'* . □ 

Let us define the modulus of continuity of PLL\^('y,a*, •) by 

Wniv)= sup {\PLLAjj,a\e)-PLLA^ij,a*,e')\}. 

\e~0'\<rj 

Lemma 2. Under SI and H1-H6, there exists a positive sequence {£k)k>i with ~^ 
when k — ?■ oo, such that for all k>l, 

."•^•■(li„s.p(»V„(i)>..))=0. 

As a consequence, the functions 9 t-^ PLL\^{'y, a* ,0) and 9 ^ Kn{9,9*) are continuous. 
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Proof. 



< Wi,n ( ^ ) + W2^n 



where 



and 



W2 



sup 

\e-e'\<i/k 



sup 

\e-e'\<i/k 



n| -/A 



exp(-/i"*'*(a;,7)) -exp(-/i"*''' {x,-i))Ax 



h"''\x,^-5,)~h<-'''\x,^-5,) 



xeW* (7)nA„ 



From SI, wc know that h°' '^{x, 7) and {x, 7) are either both finite or both infinite. 
Thus, from H6 and H3, 



< 



kj - \An\ 



15(7-2;) I da:. 



where 7^2;, for n sufficiently large denotes the configuration 7 translated by —x. Now, 
according to the ergodic theorem and by the stationarity of /i" , 



<2d 



]-)e'<E^.,0'[\9h)\]. 



From Proposition 1 and SI, h°' '^{x^-f — 6x) and h°' {x,j — Sx) are both finite when 
X e M"' (7). Thus, from H6 and H3, 



W2,n 



< 



1 



|A„| 



E 



xeTZ"' (7)nA„ 



|5((7-4)-.)|. 



Applying the ergodic theorem to the right-hand side. Proposition 2 and the stationarity 
of yu" then lead to 



W2. 



<25 - [|g(7)|c-"° ' (0,7)]<2e^J - i?„.,,.[|g(7)|] 



Therefore, for n sufficiently large 



<c5[ j: ]E^,,e'Ml)\). 



where c is a positive constant. Finally, 



limsupl Wn 



cS 



Eo.',e'i\gh)\)>ek], 
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which vanishes if one chooses, for instance, ek = 2c6{^)Ea* .e* {\g{l)\) ■ O 

4.2. Consistency of {a.n,6n) when the support parameter a* is 
unknown 

In this section, the two-step estimation procedure is considered. The consistency of this 
procedure requires some additional assumptions concerning the support parameter a. It 
must be remembered that a is a positive parameter satisfying the support hypothesis 
SI. Moreover, we assume S2-S4 below. 

52. For ah 7 e A^(M'^), a and a' mR+, OeO and A G 6(M'^), 

a<a' =^ [H^'" h) < 00 ^H^' 00]. 

53. Letting a < a* , there exists Ra > such that for every Gibbs measure /i" in 



M"*''*(^6(o*fl.)(7) = +oo)>0. 



S4. For aU 7 G A^(R'^), A e S(K'') in IR+, in 6 and a' > such that '^7) < +00, 
there exists a < a' such that -^^'^(7) < +oo. 

Assumption S2 implies that if a < a', then C A4^(R'^) and TZ'^{j) C 7^"'(7). 

As a consequence, assumption S3 appears to be a natural assumption concerning the 
support parameter. Roughly speaking, it claims that for an underestimated support of 
the energies, one could encounter some forbidden configurations with a non-negligible 
probability. 

For a given 7 G A^(R'^), S4 may be viewed as an assumption of continuity of 7^" (7) 
with respect to a. Indeed, it implies that 

]R"'(7)= U M"(7). (19) 

a<a' 

Nevertheless, in general, A4^(R'^) is not equal to Uc<<q' ■^cxjI^'')- Some consequences 
of S4 are proved in the following lemma. 

Lemma 3. Under the assumption S2 and S4, we have 

lim ^"*'«'(0^R"*-'(7 + (5o)) = 0, (20) 

lim Ai"*'^*(/i"*-''^(0,7) = +c» and 

e->-0+ , , 

(21) 

'"(0,7) < 00 andOGR" -'{j + do))=0. 
The proof of this lemma is relegated to Section 4.2.3. 
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4-2.1. Consistency of ctn 

According to S2, wc estimate the support parameter a* by the natural estimator 

a„ = inf{a > 0,H'lf{-f) < oo}. (22) 

Let us note that, from SI, ctn docs not depend on 6. Moreover, it is well defined since 
'a„ 



Proposition 3. Under S1-S3 and H1-H3, fi" -a.s., 

lim a„ = a* . 

n— ^oo 

Proof. From HI and relation (4) concerning the energies, if n < to, then 

{a,H^f{j)<^}c{a,H2f{j)<^}. 

Hence, (a„) is an increasing sequence. From (22), it is clear that d„ < a*. Therefore, 
d„ — a, where a = sup„ q;„ < a* . Let us prove that a < a* is impossible. 

Let us assume a < a* and let a' be such that a < a' < a* . Consider the average 



1 



tr„a'e / N , 1 da;. 



|A„|7a„ 

where Ra' is defined in S3. The ergodic theorem applies and for fi°' -almost every 7, 

,a* ,9' I Tja' 



lim — 

n^oo |A„ 



From S3, this last term is positive, hence for -almost every 7 and for n large 

enough, there exists x e A„ such that the energy it-/) finite. From relation 
(4), this means that for a sufficiently large tiq, '^(7) = +00. From (22), this implies 

that a' < d,i(, . But a' > a > d„f, . We have a contradiction and so, finally, -a.s., 
a = a*. □ 

4-2.2. Consistency of {a.n, On) 

The parameter 0* is estimated, as in Section 4.1, via the pseudo- likelihood procedure. 
However, a* is not known and we have to plug into definition (14) of PLL\^. The 
estimator 9n is thus defined as 

9n = argmin PLLa„ (7, d„, 6*), (23) 
see 

where is the estimator (22) and 
1 



PLLa„(7>"„,6') . 



exp(-/i"-''(a;,7))da;+ ^ /i"-^(x, 7 - 4) 



2:eR°"(7)nA„ 
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Remark 3. In the expression above, h°'"'^{x,j) is well defined only if H^^'^ij) < oo. 
From (22), this is not necessarily the case. However, one can consider a new estima- 
tor defined as an ~ an + fn, where e„ is any positive sequence which asymptotically 
vanishes. This new estimator has the same fi" -a.e. asymptotic properties as q;„ and 
/i""'^(x,7) is obviously well defined. In the sequel, we continue to denote by q;„ the es- 
timator of a* so that /i""'^(x,7) is always assumed to be well defined. In addition, if 
X S 7^"" (7), /i""'^(a;,7 — Sx) is always well defined and, moreover, it is /i" -a.e. finite 
since H'^^'^il) < oo- 

To prove the consistency of 0„, we need the hypothesis II1~H5 and the following 
modification of H6. 

H6'. One can find a real function S with S{x) — ?> when x — ?• and g G L^{n'^ ), such 
that for all a < a*, {9,9') G 9^ and 7 e X^* (R'*), if £ 7e"(7 Jq) and /i"'^(0, 7) < cx), 
then 

|/i"^'(0,7) - /i"*^'' (0,7)1 < giiMla -a*\)+ d{\e - e'\)], M"*''*-a.e. 

Let us note that in H6', we must assume that G 7?." (7 + (5o), which ensures that 
h"'^{0,j) exists when 7 G (R''). This is not a restriction in our case since we shall 
apply 116' with a = d„ (see Remark 3). Besides, from SI, the condition /i"'^(0,7) < 00 
implies that h" '^(0,7) < 00 as well. 

Theorem 2. Let fi"' eG"'-^' . Under S1-S4, H1-H5 and KG', the estimators q;„ and 
9n defined by (22) and (23), respectively, are strongly consistent, that is, /z" -a.e., 

lim (a„,4) = («*,r). (24) 

The main point of the proof of Theorem 2 is the following lemma, which is proved in 
Section 4.2.3. 

Lemma 4. Under S1~S4, H1-H4 and H6', for all OeQ, 

lim {PLLA„{l,0in,9) ~ PLLA^{j,a* ,9)) ^0, yu"*'"* -a.e. 

n— f 00 

Proof of Theorem 2. The proof follows along the same lines as the proof of Theorem 
1. Let 

K'n{9, 9* ) = PLLa„ (7, an,9)~ PLLa„ (7, «», ^*)- (25) 
Thanks to Lemmas 1 and 4, it is clear that 

lim K'n{9,9*)^Ki9,9*), 

n— >oo 

where K{-,9*) is the same function as in Lemma 1. Therefore, K'n is a proper contrast 
function. 
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Now, let the modulus of continuity of PLLa„ (7, ci„, •) be 

K(?7)= sup {|PLLA„(7,a„,0)-PLLA„(7,d„,0')l}- 
\e~e'\<r] 

The result stated in Lemma 2 for W„ is still true for W^. To prove this, it is sufficient 
to plug d„, instead of a*, into its demonstration. Hence, there exists — s- such that, 
for all fc > 1, 

(i„„„p(vn(i) 

Finally, Theorem 3.4.3 of [8] can be applied and 6'„ defined by (23) converges /i" -a.e. 
to 9. This, together with Proposition 3, completes the proof. □ 

4-2.3. Proofs of Lemmas 3 and 4- 



Proof of Lemma 3. Let us begin with the limit (20). From (19) and S2, we have 
lim M"'^'*(0^7^"*"' (7 + ^50)) 

£-i.0+ 



'■'lf]{oi 7^"*-^(7 + So)} = m"*''' U io e ^"*"^(^ + -^o)} 



^e>0 ^ ^e>0 

= M"*■^'(0^7^"*(7 + <5o))-0. 
Now, for the limit (21), 

limsupAt"''^*(/i"'"''^(0,7) = +oo and /i"*'''(0,7) < cx) and G 7^"*-'(7 + (5o)) 

<limsup/x"*''''(3Aei3(M'') s.t. e A,Hf-'-yj + Sq) = +00 

and H'^ '^(7 + (5o) < 00 and H'^ ^'^(7) < +00) 
<limsup/i"*''''(3AeS(R'^) s.t. Oe A,if^^7(o^^.(7 + <5o) = +oo 

c-i.O+ 

Thanks to a monotone class argument, 
limsupAi"'-''*(/i"'-'^''(0,7) = +oo and /i"*^^(0, 7) < cx) and e 7e"*-'(7 + <5o)) 

e-i.O+ 
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< m"* ({7 + 50 e (R")} n U n {i?^;^^(|f,i)(7 + '5o)<+oo}] 



«>OAeB(ffi<i) 

oeA 



< m"*-''* {7 + 50 e A^S: (»')} n Ui^m^f (7 + ^o) < +^} 

^ £>0 

Finally, from S4, 

limsupAi"''^*(/i"'~''''(0,7) = +(X) and /i"*'^(0,7) < 00 and G 7e"*"'(7 + (5o)) 

e-^0+ 

< /."'^^* (7 + 5o G {R") and H^l„%{j + So) = +00) = 0. 
Lemma 3 is thus proved. □ 
Proof of Lemmas 4. Let us split the difference as 

PLLa„ (7, a„, 9) - PLLa^ (7, a*, 9)^ 7^i,„ + D^^n, (26) 

where 

Di.n = j^ [ exp(-;i"-«(a;,7))-exp(-/i"*'«(a;,7))da: 

and 

^2.n=nr^ E (lK-(7)(2^)/i"''''(^>7-4)-l7^^*(^)(x)/i"*''(x,7-4)). 

In the integral of for a point x belonging to A„, there are several exclusive cases. 
The first is /),""'^(x,7) = +00 and '^(x,7) ~ +cxd. In this case, each term vanishes. 
The second is h°'"'^{x,j) < +00 and h" '^(a;,7) < +00; let us denote by Ai_„ the set of 
such x's. The last case is ft,""'®(x,7) = +00 and /i" '^(a;,7) < +00; let us denote by A2,„ 
the set of such x's. Note that because of S2, /i""'^(a;,7) < +00 and h" '^{x,j) = +00 is 
impossible since q;„ < a*. Thus, 

\D,,,\<-^ [ |e-'''-'(-.^)-e-''°*'''(-'-^)|dx+-^ / e-''°*''(-'^)dx. 

According to 113 and H6', 

/ |e-'''"'''(-.^)-e-''°*'''(^''^)|dx<e^-^ / |g(7-.)|5(|a„ - a*|)dx, 

where 7_i: denotes the configuration 7 translated by —x. If we let e > 0, then for n 
sufficiently large, from Proposition 3, 

^ / |c-'''""(-'^)-e-''°*''(-'^)|dx<.^ / |5(7-.)|dx. 
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Since g G L^{p°' ), the ergodic theorem apphes to the average m the right-hand side 
and, for n sufficiently large, 



1 



The stationarity of leads to 

/ |c-''""''(-'^)-e-''°*'''(-'^)|dx<2ee^^i?<,.,,.(|<?(7)|). (27) 



Besides, from H3 and the definition of A2,„ 

l-{ft°*.«(a;,7)<oo}-^{/j°".''(x,7)=+oo} 



If we let e > 0, then for n sufficiently large, d„ > a* — e. Hence, provided x £ ^^(7 + 
Ja;), we deduce from SI that h°'"-'^{x,j) = +00 yields ft," ~^'*(x,7) = +00. Therefore, for 
n sufficiently large. 



|A„| 



e 



A2,„ 

if 



I An I 



^ t: — r / l{/i°*.e(2;,7)<oo}n{x^K"*-'!(7+5^)} 



A„ 



|A„| 



l{;i°*.«(a;,7)<oo}n{2;G7?,°*-^(7+(5^)}n{?i°*- = .''(a:,7)=+oo} "i^- 



According to the ergodic theorem and by the stationarity of /i" , for n sufficiently 
large, 

I e-''"''(^'''Mx<2eV'''*({/i"*''(0,7)<oo}n{0^7e"*-^(7 + 5o)}) 

+ 2c^/^"'^^* ({/i"*^^(0, 7) < 00} n {0 e 7e"'-^(7 + (5o)} 
n{/."*-''^(0,7) = +oo}), 

which is less (up to 2e^) than 

/."'■^*(0^7e"*-^(7 + <5o)) 

+ //"*'^* (;i"'-^^^(o,7) = +oo|{/i"*^''(o,7) < 00} n {0 e 7e"'-^(7 + ,5o)}). 

This last term vanishes when e — >■ 0, as stated in Lemma 3. This result, together with 
(27), proves that in (26), -Di,n vanishes when n goes to +00. 
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Now let us investigate the behavior of D2.n in (26). We have 
where 



and 

^22, n 



^ 1k.„(^)(x)|/i"-^(x,7-<5.)-/i"'''(x,7-<5.)|. 



If we let e > 0, then for n sufficiently large, Q!„ > a* — e and, according to S2, 7?." ^(7) C 
71"" (7). Thus, 

^21,„<^ E (l7^=*(7)(^)-l7^"*-(7)(2^))|/^"*''(^'7-4)|. 

The application to the right-hand side of the ergodic theorem, combined with Proposition 
2 and the stationarity of /i" , leads to 

i?2i,n<2i?«^,4(l-l7^--^(7+6„)(0))|/»"*'''(0,7)|e-""*'''*(°'^)] 

<2i?„.,,4lo^^„.-.(,+,„)|/."*'''(0,7)|e-''°*''' 

Since e < e' ^ 7^"'-^'(7 + 5o) C 7^"*-<^(7 + (5o), 

SUp{lo^-R°*-«'(7+5o)J' = ^0^K°--=(7+5o) 

£'<£ 

and (20) in Lemma 3 implies that, /x" -a.e., limg_>.o 1o^-r,o*-.(^+5q) = 0. Thus, according 
to Lebesgue's dominated convergence theorem and H4, 1)21, « vanishes asymptotically. 
For D22,ni we apply H6' and for n sufficiently large, 

-D22,«<j-^ ^ lKd„(^)(x)|g((7-(5:,)_^)|. 

According to the ergodic theorem. Proposition 2 and the stationarity of ^" , we have, 
for n sufficiently large, 



^22,„<2e£;„.,e.[|5(7)|e- 



-h" (0,7)1 



which is less than 2ec^ E^* ,g-[\g{'-^)\\. Therefore, £'22,n also vanishes asymptotically and, 
as a consequence, lim„^oo -D2,n = 0. 

Returning to (26), we have proved that both Di n and 132, n vanish when n 00. □ 
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5. Examples 

In this section, we give three examples of Gibbs models which satisfy the assumptions of 
Section 4. The first is the classical model of hard spheres with pairwise step interaction 
and is hereditary. The second is the model of Delaunay tessellations with hardcore inter- 
action on the size of the cells. It is a non-hereditary model and has been recently studied 
in [5]. The third is a model where each point interacts only with its k nearest neighbors. 
A hardcore part forces these neighbors to be not too far away. These examples show that 
our assumptions are satisfied for a large class of models. 

5.1. Hard sphere model with pairwise step interaction 

Let p be a positive integer and let < ri < r2 < • ■ • < be p real numbers. For every 
a G /C (/C is a compact set in M+ \ {0}) and every 6 = {9i, . . . ,9p) in (9 is a bounded 
open set in M.P), the energy H'^' (7) is defined for every A e B{M.'^) and 7 in y\4(R'') as 

Ht^il)^ X] +ool[o,i/a](|2^-2/|) + ^il]i/a,i/Q+ri](|a:^-y|) + --- 

(28) 

Remark 4- We choose a in a compact set K. to simplify some later calculations. In 
general, for the hard sphere model, the natural hardcore parameter is the radius of the 
hard sphere (see [6]). In our case, wc take the inverse ^ so that S2 is fulfilled for a. 

Let us remark that the ranges of steps for the smooth part of the interaction depend 
on the parameter a. This seems natural since the first step occurs after the hardcore 
radius l/a (otherwise, the model could be not identifiable) and, moreover, the lengths 
of the steps, (ri, . . . , Vp), are fixed a priori. 

This model is hereditary and the existence of Gibbs measures is well known (see [12], 
for example). Let us prove that it satisfies all the assumptions needed for our estimation 
procedure. 

Proposition 4. The family of energy functions {H^^) for the hard sphere model satisfies 
assumptions S1-S4, H2~H5 and KG'. 

We first establish a lemma which is useful for checking S3. 

Lemma 5. S3, defined below, implies S3. 

S3. \/a<a*, 3Ra >0, W-f e M'£ (R"^) , 3(ai, 03, . . . , a™) GB(0,i?a)™ (m>l),3e> 
0, Vxi e i3(ai,e),Vx2 G i?(a2,e), ■ • ■ ,Vxm G B{am,s), 

^b(q,b.^)^^b{o,r^Y 1-4„) < +00 
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and 



H 



B{0,Ra 



)(7S(0,. 



+00. 



Proof. Suppose that (H'^' ) satisfies S3. Let a < a* and /x" be a Gibbs measure in 



(d7) 



{«2(0,H„,(7..(0,«„,e+7')= + Oo} >^*(75(0,K„)c) 

From S3; for almost every 7 in (K''), 

/■ g--f^B(o!'n„)(TB(o,R„)'=+')'') 

Therefore, /x" {^b'(o r. = +00) > 0. The lemma is thus proved 



-^B(o,fl„)(d7')A^"''' (d7)- 



7rs(o,i?„)(d7') > 0. 



□ 



Proof of Proposition 4. SI and S2 are obviously satisfied. The choice of ^ instead of 
a for the hardcore radius parameter is crucial here. 

Choosing i?„ = |, r = - i), m = 2, oi = (0,0) and 03 = + 2r,0), we check 
that S3 in Lemma 5 is satisfied, which implies S3. 

S4 is satisfied because the number of terms in (28) is finite. So, given 7 such that 

iJ^ '^(7) < +00 for some a' in /C, it is easy to find a < a' such that H'^'^['y) remains 
finite. 

Moreover, as the energy is clearly invariant by translation, H2 is obvious. Let us note 
that for every 7 e M'^{R'^), 



^^'(0,7) 

= +OOl[o,i/a] i\x\) + 6'll]l/a4/a+n] (kl) + ' ' ' + dp^]l/ a+r^^^ ,1/ a+r^]i\x\) ■ 

xe7 



(29) 



It is easy to see that /i"'^(0,7) is either equal to infinity or uniformly bounded. Conse- 
quently, assumptions H3 and H4 are satisfied. 

H5 is deduced from a classical identification property of the exponential family models 
(see, for example, assumption [ident] on page 244 of [2]). 

It remains to prove 116'. For every 7 in A^J^ {^'^), we define 



d{'y) = min 



r2 



1 
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d(7) is the minimal distance, over every x in 7, between \x\ and a discontinuity point of 
the pairwise potential in (29). Since /i" is locally absolutely continuous with respect 
to the Poisson process, the random variable d is fi" -a.s. positive. Consequently, there 
exists a positive function ip on R+ \ {0} such that ip is decreasing, lun^^Q+ ip(x) = +00 
and ip{d{-)) is in L^ifJ." )■ It suffices to choose, for instance, = where F is the 

repartition function of d(-) for the probability fj," . 

From the expression (29) of h"' {0,j), there exists K >0 such that, for every a < a* 

and 7 in A4'^ (W^) which satisfy S TZ" (7 + So) and /i"'^(0, 7) < +00, we have, for every 

{9,0') in 92, 



|/i"''(0,7)-/i"*'' (0,7)1 <i^l[d(^),+oo[ 



1 1 



< K- 



(p{l/a — 1/a*) 
< i^max((p(d(7)),l) 



~K\0-e'\ 

K\e-e'\ 
1 

p{K'\a — a* I) 



+ \e-e'\ 



where K' is a positive constant. This inequality implies H6'. Proposition 4 is thus 
proved. □ 



5.2. Hardcore Delaunay tessellations 

In this section, we present the model of Delaunay tessellations studied in [5] . A hardcore 
interaction forces the cells of the tessellation to be not too small and not too large. Let 
us state some definitions concerning Delaunay tessellations. We assume that d = 2. Let 
7 be in A4{M.'^). A subset T of three points in 7 is a Delaunay triangle if the open 
circumscribed ball B{T) of T does not contain any point of 7. Del{'~f) denotes the set of 
Delaunay triangles of 7. 

Let r > be a positive real number. For every a € /C (/C is a compact set in M+ \ {0}) 
and every 6* e 6 (9 is a bounded open set in E), we define the energy if^'^(7), for every 
A in B(R'^) and 7 in M{R'^), as 

E ^(^) (30) 

TeDei(7) 
B(T)nA5<i0 

with 

( +00, if i?(T) > a, 

ip{T) = l +00, if l{T)<r, (31) 

[ ePer{T), otherwise, 

where R{T) is the radius of B{T), 1{T) is the smallest length of the sides of T and Per{T) 
is the perimeter of T. 
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Remark 5. In this example, we choose to estimate the hardcore parameter a (as the 
maximum radius R{T)). It would also be possible to estimate the parameter r (as the 
minimal length l{t)). Moreover, we believe that it is possible to estimate both of these 
parameters simultaneously, but this situation is not really dealt with in this paper and 
some further assumptions and proofs would be needed to deal with this case. 

We have chosen the smooth part of the energy to be equal to 9Per{T) because it seems 
interesting to us. However, other smooth functions depending on the volume, the angles, 
etcetera are valuable and they could be managed in the same fashion. 

Proposition 5. The family of energy functions {H^'^) for the hardcore Delaunay tes- 
sellations model satisfies assumptions S1-S4, II2-II5 and H6' if 

"iaeK a>r. (32) 

Proof. SI, S2 and S4 are obviously satisfied, as in the example of Section 5.1. To prove 
S3, we use Lemma 5 and must prove that S3 is satisfied. 

Lemma 6. S3 is satisfied. 

Proof. Assumption (32) is crucial here. Let be in and a < a*. Let T ~ {01,02,02} 
be an equilateral triangle in M**, where the radius R{T) of the circumscribed ball B{T) 
satisfies a < R{f) < a*. Moreover, let Ra > he such that B{f) C B{0,Ra - r - a*). 
Now, for every 7 in A4^ (R''), we must add some points 04, 05, ... , o,„ to the points oi, 
02, 03 to obtain the expected configuration needed in assumption S3. 

We construct these points 04, 05,..., o™ recursively. This construction is inspired by 
the one given in [5], page 148. Let us start with {04,05, ... = 7B(o.R„)\B(o,flc-a*) 

and put 7 = ^B{o,Ra)'' + + + + ^ h Sa^, ■ Wc now test whether 7 satisfies 

the assumption 

V.TeB(0,i?„) j{Bix,a*))>0, (A) 

where B{x,q*) denotes the closed ball centered at x with radius a* . If this is the case, 
we stop the recursive process, put m = m' and oi, 02, ... , o^ are obtained. If (A) is false, 
we choose an arbitrary point Om'+i in B(0,Ra) such that 7(i?(o,„'+i,a*)) = and we 
put 7 := 7 + Sa^,^-^- We then test (A) for this new 7. If (A) is true, we stop here and 
TO = to' + 1; otherwise, we choose another point a„i'+2 as above and put 7 := 7 + Sa^,^.^ ■ 
We continue like this until the process stops, which is always the case since B{0, R^) is a 
bounded set and, from (32), the minimal distance between two points in 7 is larger than 
r. Finally, it is relatively easy to see that 

^b(o,R^)^'^B(0,R^Y + <5ai + h 5a,J) < +00 

and 

^B(0,_R„)(TB(0,fl»)= +Sai-\ +'5a„) = +0O. 
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Both of these properties remain true for small local perturbations of ai, a2, . . . ,a„i. So, 
for £ sufficiently small, S3 is proved. □ 

H2~H5 and H6' still have to be proven. The energy is clearly invariant under transla- 
tion, so H2 is obviously satisfied. Moreover, for every 7 e M'^{M.'^), 

/^"''(0,7)= E ^(^)- E ^(^)' (33) 

Te_Dei(7+<5o) TeDel(t) 
06T OeB(T) 

SO ft,"-^(0,7) is either equal to infinity or uniformly bounded (see [5], page 137). Conse- 
quently, assumptions H3 and H4 are satisfied. Contrary to the example of Section 5.1, 
assumption H6' is obviously satisfied here. H5 still has to be proven. With a proof similar 
to the one given in Lemma 6 above, we show that V9 eQ\{9*}, 3Ra* > 0, V7 G A1^' (R'') , 
3{ai,a2, ■ ■ . G 5(0, i?Q,.)™, 3e > 0, Vxi e i?(ai,£),V.T2 G 5(02, e), . . . ,Va;,„ G i3(a,„,£) 

/i"*''* (0,7s(o,fl„.)^ + d^.i + ■ • ■ + 4,J ^ /i"'''(0,7B(o,i?„.)'= + S.J. (34) 

With an argument similar to that used for Lemma 5, it can be proven that (34) implies 
H5. Proposition 5 is thus proved. □ 



5.3. A forced clustering fc-nearest-neighbors model 



In this section, wc introduce a model where any point x interacts with only its k nearest 
neighbors (sec [1], for example). We propose including a hardcore interaction which forces 
these k nearest neighbors to be at a distance less than a from x. Let us describe the 
interaction precisely. 

Let fc > 1 be a fixed integer and let ip he a bounded function from M+ to R which is 
non-null in a neighborhood of 0. 

For every a G /C (/C is a compact set in R+ \ {0}) and every G (9 is a bounded 
open set in R), we define the energy iJ^'^(7), for every A G B(R'^) and every 7 in M{W'-), 
as 



(7)= E +oot[o^k]{l{B{x,a)))+ M\x-y\) 



(35) 



where A" denotes the set \J^^f^B{x,a), N^[x,^) is the set containing the k nearest 
neighbors of x in 7 and B{x^a) is the closed ball centered at x with radius a. 

The hardcore interaction implies that the locally finite energy configurations do not 
contain clusters with less than k+\ points. This model therefore naturally forces clusters 
and is a non-hereditary model. 

In the literature, no proof for the existence of Gibbs measures for these kinds of inter- 
actions is available, but we claim that the proof is essentially the same as that given for 
Example 2 in [5]. 



Non-hereditary point processes 



1395 



Proposition 6. The family of energy functions {H'^' ) for the forced clustering k- 
nearest-neighbors model satisfies assumptions S1-S4, H2-H5 and H6'. 

Proof. SI, S2 and S4 are satisfied as in both the examples of Sections 5.1 and 5.2. S3 is 
satisfied by checking S3 in Lemma 5, which can be easily done. For the other assumptions, 
note that for every 7 S M'^{M.'^), 

/i"'«(0,7) = +ool[o,fc](7(S(0,a)))+ ^ 9ipi\x\) 

(36) 

+ MN)-Mlx'(^,7)l), 

where x'^i^il) is the farthest neighbor of x in Af^i^x,^). 

In (36), the number of terms in the first sum is obviously equal to k. In the second sum, 
the number of terms is uniformly bounded (for instance, by bk when d — 2). Consequently, 
as in both the examples of Sections 5.1 and 5.2, /i"'^(0,7) is either equal to infinity or 
uniformly bounded. Therefore, assumptions H3, 114 and H6' are satisfied and 115 is proved 
by an argument similar to that used in the example of Section 5.2. □ 
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